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Lattice QCD using fermions whose Dirac operator obeys the Ginsparg- Wilson relation, is perhaps 
the best known formulation of QCD with a finite cutoff. It reproduces all the low energy QCD 
phenomenology associated with chiral symmetry at finite lattice spacings. In particular it explains 
the origin of massless pions due to spontaneous chiral symmetry breaking and leads to new ways 
to approach the (7(1) problem on the lattice. Here we show these results in the path integral 
formulation and derive for the first time in lattice QCD a known formal continuum relation between 
the chiral condensate and the topological susceptibility. This relation leads to predictions for the 
critical behavior of the topological susceptibility near the phase transition and can now be checked 
in Monte-Carlo simulations even at finite lattice spacings. 
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I. INTRODUCTION 

Formulating QCD non-perturbatively is important 
for understanding the low energy phenomenology of 
hadronic physics. The only such formulation that is 
known at present is through the lattice regularization. 
However, until now fermions have been a major obstacle. 
The main reason has been that chiral symmetry is eas- 
ily broken on the lattice. Since two of the main features 
of low energy QCD, namely the spontaneous breaking 
of chiral symmetry and the anomaly, are intimately con- 
nected with chiral symmetry, one hardly uses the lattice 
discretization to discuss the low energy phenomenology 
of QCD. Even if such a discussion was carried out, one 
usually needs to invoke continuum limits before the for- 
mulation can actually be argued to reproduce all the low 
energy properties of QCD. 

At finite lattice spacings it is usually unclear as to how 
chiral symmetry is realized and what properties of QCD 
are lost. For example with Wilson fermions [|l|, since chi- 
ral symmetry is explicitly broken, there is no reason to 
expect massless pions unless one tunes a mass parameter 
to a critical point. Further, the critical point appears to 
describe spontaneous breaking of parity and flavor sym- 
metry 1^. Thus the soft pion theorems of QCD will not 
be reproduced at finite lattice spacings. If instead one 
starts with staggered fermions ||^ one can only formu- 
late four flavors of QCD in the Lagrangian formulation. 
In this case a U{1) subgroup of the flavor non-singlet 
chiral symmetry is exact on the lattice. However, the fla- 
vor symmetry gets broken at flnite lattice spacings which 
again is expected to be restored only in the continuum 
limit. 

The anomalous breaking of the flavor singlet chiral 
symmetry is expected to give the rj' particle its mass. 
With both Wilson and staggered fermions, discussions in- 
volving the anomaly have new complications. Since the 
lattice regularization breaks the anomalous chiral sym- 
metry along with other chiral symmetries, the effects of 
the anomaly cannot be easily isolated. The simple con- 



tinuum discussion in terms of the zero modes of the Dirac 
operator and topology of the gauge flelds breaks down. 
All this indicates that both Wilson and staggered fermion 
formulations are not ideal for a description of low energy 
QCD phenomenology at finite lattice spacings. 

A long time ago Ginsparg and Wilson Q showed that 
if the Dirac operator, D, obeyed the relation 



■D75 + I5D = aDj^D, 



(1) 



where a is the lattice spacing, then the theory would 
have a remnant chiral symmetry up to contact terms. An 
year ago Hasenfratz noticed that the fixed point action of 
QCD obeyed the Ginsparg- Wilson relation Q and hence 
a series of interesting results [^,0 about renormalization 
and lattice index theorems followed. The new massless 
fermion formulations motivated by the overlap Q also 
obeyed the Ginsparg Wilson relation. It seemed that 
there was more to the relation than was appreciated in 
the past. Recently Liischer showed that there is an exact 
chiral symmetry on the lattice |^ due to the Ginsparg- 
Wilson relation, since the Nielsen-Ninomiya theorem |1C| ] 
is not applicable. As a result, for the first time, we can 
write down a formulation of QCD with the right inter- 
nal symmetries of the theory at finite lattice spacings. 
Thus in this formulation the role of the continuum limit is 
merely to renormalize some multiplicative constants and 
suppress finite lattice artifacts which do not change the 
internal symmetry structure of the theory. This leads to 
the most elegant non-perturbative formulation of QCD. 

In this paper we review the main features of the above 
observations. We use the exact chiral symmetry that 
arises due to the Ginsparg- Wilson relation to derive the 
essential physics of low energy QCD in the path integral 
formulation. For the first time we are able to discuss 
the physics of chiral symmetry breaking relevant to QCD 
with any number of flavors, the associated pion physics 
and the essential ingredients of the U{1) problem, all di- 
rectly at finite lattice spacings without much effort. In 
particular we derive a chiral condensate of the theory 
which reflects the breaking of chiral symmetry. For two 



1 



or more flavors we reproduce the Goldstone's theorem if 
the chiral symmetry is spontaneously broken. Further, 
since the Ginsparg- Wilson relation allows for the break- 
ing of the anomalous chiral symmetry through exact zero 
modes of the Dirac operator like in the continuum, we 
can easily discuss a solution to the U{1) problem. By 
introducing an explicit quark mass we can derive a well 
known continuum relation, between the chiral conden- 
sate and the topological susceptibility, in lattice QCD. 
The problems associated with the quenched approxima- 
tion can be traced to the violation of this relation. Near 
the chiral phase transition at finite temperatures, this 
relation also leads to predictions for the critical behav- 
ior of the topological susceptibility which have not been 
appreciated before. These predictions can be verified in 
Monte-Carlo simulations. 

Most of the results that we derive already appear in 
However, we have included them for the coher- 
ent presentation of our work, and as far as we know this 
is the first time the results have been derived in the path 
integral formulation based on the symmetry discovered 
by Liischer |^ . This makes the re-derivations sufficiently 
different and perhaps more transparent. Our presenta- 
tion of the C/(l) problem and its consequences are new. 
We have tried to keep the derivations simple and easily 
accessible to readers familiar with field theory. 

II. CHIRAL SYMMETRY 

Let us begin by reviewing the exact lattice chiral sym- 
metry, discovered by Liischer that arises due to the 
Ginsparg- Wilson relation. The partition function of lat- 
tice QCD is given by 

Z = J [d[/]e"^^^f^l y[d^][d*]e(-^^t*'*l\ (2) 

where SgIC^] is some lattice gauge action and is the 
bare gauge coupling. The fermionic action is given by 

Sf[^, = -a^ J2 ^{x)D^,y^{y), (3) 

a:, J/ 

where we assume that D obeys eq.(|]). We have sup- 
pressed color and fiavor indices for convenience. The 
action is invariant under the infinitesimal chiral trans- 
formations of the form 

= £75(1 - a^D)^, (5* = e*(l - a^D)j5. (4) 

However, the measure is not invariant under these trans- 
formations. In fact if = [del)] 

S{d[(l)]) = eatT{-f5D)[d(l)] = -~2eNf{n+ - (5) 

where Nf is the number of quark flavors |9j . This means 
that the flavor-singlet chiral transformation suffers an 



anomaly on the lattice like in the continuum. Further 
like in the continuum the anomaly arises only due to a 
non-zero index of D which is given by (n+ — n_ ) , where 
n± are the number of zero modes of D which are also 
simultaneous eigenstates of 75 with eigenvalue ±1. We 
will explain the origin of cq.(|5|) later. The flavored chi- 
ral transformations are obtained in the usual way by in- 
cluding a traceless flavor matrix in the transformations 
of eq.(p. The flavored transformations are exact sym- 
metries of the theory. Thus chiral symmetry and the 
anomaly emerge exactly like in the continuum, however 
all the above arguments are being made on a finite lat- 
tice. 

A natural candidate for the chiral condensate emerges 
like in the continuum. Starting from the definition of the 
(unnormalized) expectation value, (^'754') and perform- 
ing a change of integration variables by an infinitesimal 
chiral transformation, it is easy to show that 

(*75*) = (*75*) + 2e(*(l - ^aD)^'). (6) 

This means that if the theory has an exact chiral sym- 
metry then 

(^(l - ^aD)^) = 0. (7) 

When the theory contains more than one flavor, the 
above relation strictly holds at flnite volumes. This is 
because there exist flavored (chiral) symmetries that are 
not broken by any anomalies. The only way the chi- 
ral condensate could be non-zero, is due to spontaneous 
breaking of chiral symmetry which is a subtle feature of 
the thermodynamic limit. To see this one has to add 
a small mass term to the Dirac operator and break the 
symmetry explicitly. If the symmetry is spontaneously 
broken then the condensate would be non-zero when the 
thermodynamic limit is taken before the explicit symme- 
try breaking mass is set to zero. In the case of a single 
flavor, the only chiral symmetry of the action is broken 
by the measure and (f'(l — ^aD)"^) ^ even in flnite 
volumes. 



III. PIONS AND THE U(l) PROBLEM 

In the presence of a mass term we can derive the mod- 
ified relation 

{1 + am)[{D + m)-f5 + -f5{D + m)] 

= m{2 + am)75 + a{D + ni)j5{D + m), (8) 

using the Ginsparg- Wilson relation. To see the origin 
of massless pions when the chiral symmetry is sponta- 
neously broken, let us look at the zero momentum pion 
correlation function Gab{p = 0), which is given by 

Gab = ^5](§zT''75^'(a;)^iA5*(y)), (9) 

x,y 
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where f2a'* is the space-time volume. r°'^ are the Nj — 1 
hermitian, traceless flavor matrices. Doing the Grass- 
mann integral the above correlator can be rewritten as 
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ria^ \ \ D + m D + m 



l5-p-, 75 



(10) 



where we have absorbed the fermion determinant in the 
path integral measure along with the Boltzmann weight 
of the gauge fields used to calculate the expectation value. 
By multiplying eq.(||) with {D + m)~^ from left and {D + 
m)~^75 from the right we get 



tr 



1 



-75- 



D + m""D 
2(1 + am) 



-tr 



m 
1 



-75 



1 



m{2 + am) \D + jn 2(1 + am) 



(11) 



which can be used to rewrite Gab as an expectation value 
of the chiral condensate. In the chiral limit this leads to 



lim Gab = Sab^im - laD]^). 
m—>o ma^ I 



(12) 



If chiral symmetry is spontaneously broken the right 
hand side is singular in the chiral limit. This singularity 
shows that oc m, since in the chiral limit we expect 
Gab oc (mj)^^. This is the Goldstone's theorem. 

It is straight forward to obtain a solution to the U{\) 
problem in the above language. Consider the zero mo- 
mentum flavor-singlet pion correlator 



" ^Xl^* *'^5*(^)* *^5*(2/))- 



(13) 



The Grassmann integral now yields two terms 
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(14) 



where again we have absorbed the fermion determinant 
into the Boltzmann weight used to calculate the expecta- 
tion value. The first term is familiar from the non-singlet 
pion correlator. The second term is non-zero only in con- 
figurations that have a non-zero index of the lattice Dirac 
operator, exactly like in the continuum. This is because 
we can use eq.(|8|) and eq.(||) to show that 



tr 



1 



+ m 



75 



27V 



/ 



myl + am 



■(n_|_ — nJ) 



(15) 



is even on a finite lattice. Thus we get 



Goo = 4^i-^(*[l 



m(l -|- am/2) 



2(1 -I- am) 



{D + m)]^) 



m2(2H-am)2 n 



-((n+-n_)2). 



(16) 



This relation can be used to obtain new insight on the 
U{1) problem. A solution to the U{\) problem requires 
Goo to remain non-singular in the chiral limit. It is clear 
that the first term on the right hand side is singular in 
the chiral limit, if (*[1 - \aD]^) 7^ 0. This is expected 
to be true in finite volumes for one flavor QCD due to the 
anomaly, and infinite volumes for two or more flavor QCD 
due to spontaneous chiral symmetry breaking. In both 
these cases, if Gqo must remain non-singular, the singular 
piece in the first term must cancel with a corresponding 
piece in the second term. This requires 



((n+-n_)2) 



(*[1 - -aD]*) 



(17) 



at the leading order in the chiral limit. Thus we see 
that the solution to the U{\) problem requires a connec- 
tion between the lattice topological susceptibility defined 
using the index of the lattice Dirac operator, and the chi- 
ral condensate in the chiral limit. Such a connection is 
known in the literature in the context of effective models 
|]l2| and the large Nc limit jl^. However, now the same 
result emerges directly from lattice QCD. 

At finite volumes in the one fiavor case, eg. ([l7|) is easy 
to understand since only sectors with unit topological 
charge contribute to both sides of the equation. How- 
ever, if one takes the infinite volume limit before tak- 
ing the chiral limit, the relation is somewhat surprising. 
In this limit, the chiral condensate is expected to get 
contributions from a non-zero density of small eigenval- 
ues of the Dirac operator. This is the content of the 
Banks-Casher formula |]l^ and will be derived later in 
the present context. The exact zero modes of the Dirac 
operator on the other hand are expected to be suppressed 
with the volume. Thus it is unclear why the topological 
susceptibility, which is related to exact zero modes on 
the lattice, is also related to the density of approximate 
zero modes as eq.(|l^) requires. Another interesting con- 
sequence of the equation is that the topological suscep- 
tibility goes linearly to zero with the mass, independent 
of the number of flavors, if chiral symmetry is sponta- 
neously broken. This appears to be a consequence of the 
fact that the thermodynamic limit and the chiral limit 
do not commute fl^ . 

The relation described by eq.(^^ is valid only in the 
full theory of QCD where fermions can be created and 
destroyed in the vacuum. In the quenched limit the left 
hand side is independent of the mass and number of 
fermion flavors, which only enter through the fermion 
determinant. Thus in the quenched theory the flavor- 
singlet pion correlator will show divergences in the chiral 
limit, since the singular pieces between the two terms 
in eq.(p^) do not cancel. The form of these singulari- 
ties have been predicted using chiral perturbation theory 
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|g5[ and can now be tested even at finite lattice spacings. 
Further, since exact zero modes of the Dirac operator are 
necessary for the anomaly, they are physical and occur 
even at finite lattice spacings. These produce new sin- 
gularities in the quenched theory at finite volumes. In 
particular the chiral condensate is expected to diverge in 
the chiral limit due to the zero modes. 



IV. CRITICAL BEHAVIOR OF ZERO MODES 

We can extend the above analysis to finite tempera- 
tures. It is clear that eq.(^^ will be valid at non-zero 
temperatures in the phase where chiral symmetry re- 
mains broken and if the fiavor singlet pseudo-scalar corre- 
lator is free of infra-red divergences. It is usually believed 
that for two massless flavors, infra-red divergences arise 
at the chiral transition T^, only due to diverging corre- 
lation lengths in the pseudo-scalar isovector(the pions) 
channel and the scalar isocalar channel(the sigma). This 
leads to the 0{4) critical behavior which is second order 
[ p^ . If this is true, the fiavor singlet correlator given in 
eq.(|l^), must be free of infra-red singularities even at Tc- 
The chiral condensate is the order parameter of the phase 
transition and at Tc has a critical behavior given by 



(^-[1 - -aZ?]*) - 



(18) 



where 5 = 4.82 ± 0.05 |17| is an 0(4) critical exponent. 
The cancellation of the singular pieces in the right hand 
side of eq.(|6h leads to 



(19) 



at Tc- This critical behavior of the topological suscepti- 
bility has been ignored in the literature and can form a 
useful consistency check of 0(4) critical behavior. 

For T > Tc, the topological susceptibility need not 
be related to the chiral condensate since in this phase 
both terms in the eq.(|l^) are finite in the chiral limit. 
However, if we assume that the thermodynamic and the 
chiral limits commute in the chirally symmetric phase, 
then for two or more flavors we can also predict 



T < Te 

T > T, 



(20) 



The prediction for T > Tc is easily derivable from gen- 
eral considerations |^ . This behavior is consistent with 
what we know about the phase transition. For example, 
with three or more flavors the transition is expected to 
be first order. This is refiected in the dramatic change 
in the small mass dependence of the topological suscep- 
tibility. The change in the two fiavor theory is smaller 
and matches with the prediction of eq.(19) for a second 
order transition. 

It is useful to look for the above critical behavior of 
the topological susceptibility in any study of the QCD 



phase transition, since it refiects the interesting interplay 
between the chiral symmetry and the U{1) anomaly. It 
has been appreciated that this interplay is important in 
determining the order of the phase transition for two light 
flavors ||l^. In particular if the anomalous U{1) is small 
for some dynamical reason, the fluctuations from new 
degrees of freedom can make the transition flrst order. 
On the other hand if the symmetry is broken by a large 
amount then the phase transition is expected to be in the 
0(4) universality class. 

The fermionic definition of the topological suscepti- 
bility, used in the above analysis can in principle be 
measured on the lattice since it involves calculating zero 
modes of the Dirac operator in a realistic simulation. In 
the staggered fermion formulation, which is most suited 
for studying the chiral phase transition, zero modes are 
shifted due to finite lattice spacing effects j2^]. Such ef- 
fects distort the critical behavior of the zero modes, and 
one needs to invoke the continuum limit to see the effects 
of the anomalous symmetry emerge at the QCD phase 
transition [ pT|JT^ . However, approaching the continuum 
limit is quite difficult. On the contrary, there is evidence 
p2| that with Ginsparg- Wilson fermions the continuum 
limit is not necessary to see the zero modes that are re- 
sponsible for the anomalous symmetry breaking. Thus 
the above critical behavior should emerge even at finite 
lattice spacings. 



V. DISCUSSION AND CONCLUSIONS 

All the above results have been derived without spec- 
ifying an explicit realization of the the Dirac operator 
that obeys the Ginsparg- Wilson relation. However, it is 
important that the Dirac operator does not suffer from 
doubling, since it is possible to find Dirac operators which 
obey the Ginsparg- Wilson relation but do not solve the 
fermion doubling problem. In such cases, the real dynam- 
ics of QCD will not be reproduced. The chiral symmetry 
discussed could be fiavored as in staggered fermions. The 
two known classes of solutions based on the fixed point 
actions |^J2^ and the overlap formulation ||] are expected 
to solve the doubling problem completely. The proposal 
by Neuberger has a simple structure based on the four 
dimensional Wilson like Dirac operator Dw , and is given 
by aD = 1 + Dw{dI^Dw)~^^'^ ■ The important point 
is that Dw must have an appropriately tuned negative 
mass term to be in the right phase to produce QCD. 
Whether such a phase exists is currently being explored. 
A further effect of the Wilson mass is to renormalize the 
theory and control leading lattice artifacts. 

The Neuberger Dirac operator is also closely related 
to the Shamir's variant of the domain wall fermions [ p5[ , 
when the distance between the walls is taken to infln- 
ity. It is likely that all the nice properties discussed 
above will emerge with small violations when the walls 
are separated by a large but flnite separation. This strat- 



4 



egy can be used for practical simulations and still retain 
the interesting properties of Ginsparg- Wilson fermions. 
Preliminary tests for QCD has proved that domain wall 
fermions preserve chiral symmetry properties very well 
p6t . Presently, these Domain wall fermions are being 
used to study the QCD phase transition |27|. Prelimi- 
nary tests seem to indicate the presence of zero modes 
that break the anomalous symmetry even above 

the phase transition. This suggests that it should be 
possible to extend the above discussion of the critical be- 
havior of the topological charge to such a study. This 
would give a more quantitative grasp of the two flavor 
QCD phase transition than has been possible before. 

The above solutions to the Ginsparg- Wilson relation 
obey an extra hermiticity property, namely Z?^ — 75-D75. 
Based on this it is possible to write a formula for the 
chiral condensate in terms of the eigenvalues of the Dirac 
operator. It is easy to show that aD — 1 is unitary whose 
eigenvalues can be represented as — exp(z0), — tt < < 
TT. The exact zero modes of D correspond to eigenvalues 
at 61 = 0. Thus 



d0 p{e) 



1 + cos{9) + i sm{e) 
2(1 - cos(6') - i sm{0) + am) ' 



(21) 



where p{9) is the normalized density of eigenvalues 
/ d9 p{e) = ANjNc. Due to the hermiticity of D the 
eigenvalues come in complex conjugate pairs, so that p{0) 
is a symmetric function. In the m ^ limit it is easy to 
show that 



(§[1 - \aD]^) = Jj^p(O). 



(22) 



This is the Banks Casher formula In the present 

context it is easy to see that tr(75_D) gets contribution 
only from the subspace of eigenvectors of D with eigen- 
values 2/a. However this contribution is also related to 
(n+ — n_) due to the fact that tr(75) = 0. This is one 
way to understand the origin of eq. (H) . 

Results from simulations using the Dirac operator sug- 
gested by Neuberger in the Schwinger model |Q and the 
domain wall fermions in quenched QCD [ p9| give clear ev- 
idence for divergences in the quenched chiral condensate 
due to exact zero modes on finite lattices. Such singular- 
ities are expected due to the anomaly as discussed above, 
but have not been seen in realistic simulations involving 
staggered fermions. This perhaps reflects the fact that 
the breaking of chiral symmetry in the quenched theory, 
is represented well by Ginsparg- Wilson fermions even at 
finite lattice spacings. Thus, quenched lattice QCD with 
Ginsparg- Wilson fermions may give a more definitive an- 
swer to the question about the reliability of the quenched 
approximation. 

In conclusion, we have shown how Ginsparg- Wilson 
fermions reproduce the low energy physics of QCD at 
finite lattice spacings. The exact chiral symmetry of 



the action can in principle be used to derive a host of 
soft pion theorems even at finite lattice spacings. We 
have concentrated on the Goldstone's theorem and ob- 
tained interesting consequences from a solution to the 
?7(1) problem. We have also suggested how the critical 
behavior of the topological susceptibility can confirm the 
second order nature of the QCD phase transition with 
two massless quarks. If simulations involving Ginsparg- 
Wilson fermions are feasible it is very likely that the re- 
sults would reflect physical reality more easily than the 
conventional fermion formulations. 
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